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Qh. Measurements are presented of the polarisation of W + W~ boson pairs produced in e + e~ collisions, 
and of CP-violating WWZ and WW7 trilinear gauge couplings. The data were recorded by the OPAL 
experiment at LEP during 1998, where a total integrated luminosity of 183 pb _1 was obtained at a 



centre-of-mass energy of 189 GeV. The measurements are performed through a spin density matrix 
analysis of the W boson decay products. The fraction of W bosons produced with longitudinal 
polarisation was found to be OL/Vtotai = (21.0 ± 3.3 ± 1.6)% where the first error is statistical and 
the second systematic. The joint W boson pair production fractions were found to be orr/^totai = 
(78.1 ± 9.0 ± 3.2)%, <7 L i> total = (20.1 ± 7.2 ± 1.8)% and <7 T i>total = (1.8 ± 14.7 ± 3.8)%._ In 
the CP-violating trilinear gauge coupling sector we find k z = — 0.20l 007 , g\ = — 0.02]to'33 and A z = 
— 0.18l ' 16 , where errors include both statistical and systematic uncertainties. In each case the coupling 
is determined with all other couplings set to their Standard Model values except those related to the 
measured coupling via S77(2)l x C(1)y symmetry. These results are consistent with Standard Model 
expectations. 
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1 Introduction 



We report measurements of the properties of W pair production in e + e~ collisions using data recorded 
by the OPAL detector at LEP at a centre-of-mass energy of 189 GeV with a total integrated luminosity 
of 183 pb _1 . We perform a spin density matrix (SDM) analysis Jl],§] of the production and decay 
properties of the W bosons using the semi-leptonic WW — ► qqeu, WW — ► qq/xz/ and WW — ► qqru 
final states. Using suitable summations of SDM elements we present measurements of the inclusive 
production cross-sections for each of the transverse and longitudinal polarisation states of the W. The 
SDM elements are also sensitive to triple gauge couplings. We present measurements of CP-violating 
couplings involving the W 1 * 1 , Z° and photon. 

The doubly resonant e + e~ — > W + W~ production process proceeds via s-channel Z° or photon 
exchange, or via t-channel neutrino exchange, collectively known as the CC03 diagrams. The s-channel 
processes contain the WWZ and WW7 triple gauge boson vertices. The most general Lorentz invariant 
Lagrangian describing these vertices ]3|,||| contains 14 independent couplings. In order to facilitate 
measurements the number of parameters is often reduced to three by assuming SU{2)\ J x U{l)y gauge 
invariance and charge conjugation (C) and parity (P) invariance. The resulting independent couplings 
are conventionally taken as: Ak 7 , Agf and A ||. Measurements of these couplings using data recorded 
at LEP2 |]-10| and the Tevatron [11] have been reported elsewhere. 
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If C, P and CP-invariance are not assumed then several additional couplings may be present. The 
CP-violating ones can be taken as ky and Ay which violate P and conserve C, and g\ which violates 
C but conserves P ( V = Z° or 7) |J. A further parameter, glf , violates both C and P but conserves 
CP so is not considered here. SU (2)i, X U(l)y symmetry imposes the following relations [12-14] which 
are assumed in this analysis: 



(1) 



All of the triple gauge boson couplings (TGCs) can in principle be measured through observations 
of the W production angular distribution, and distributions of the W decay products. One way of 
realising such an analysis is through the spin density matrix (SDM) approach. In this approach the 
individual contributions to the W production angular distribution arising from each of the different 
possible helicity states of the W bosons can in principle be determined exclusively. These exclusive 
SDM distributions exhibit different behaviour with respect to each of the TGCs. The TGCs can 
therefore be determined from the SDM elements in a second step. In W boson pair production, 
the SDM analysis is particularly suited to the extraction of the CP-violating couplings. Indeed, 
the imaginary parts of the off-diagonal elements of the SDM are completely insensitive to the CP- 
conserving couplings, and will only deviate from their Standard Model predictions in the presence of 
CP-violation at the triple gauge boson vertex. 

Before proceeding to the TGC measurements we investigate the exclusive production cross-sections 
for each of the transverse and longitudinal polarisations of W bosons. These can be made in the context 
of the SDM analysis by suitable summations of SDM elements which are described in detail in the 
following section. The study of the longitudinal cross-section is particularly interesting as this degree 
of freedom of the W only arises in the Standard Model through the electro-weak symmetry breaking 
mechanism. Previous measurements of the proportion of longitudinally polarised W bosons produced 
at LEP2 have been been reported by the OPAL || and L3 (l5| experiments. 

In addition we present limits on the CP-violating parameters. Previous limits on the CP-violating 



TGCs have been obtained at the DELPHI experiment using data at 161 GeV and 172 GeV ||. Limits 
on the CP-violating WW7 couplings, k 7 and A 7 , have been obtained at the DO experiment [jl6| from 
the process pp — > luj + X. Measurements of the neutron electric dipole moment show that the 
electromagnetic interaction is CP-conserving to very high accuracy [17]. 



2 The Spin Density Matrix 

The two-particle joint spin density matrix (SDM) Q describes the contribution of each of the helicity 
states of the W bosons to the W + W~ production cross-section. 

The amplitude for the production of a W~ with helicity r_ (= —1,0 or 1) and a W + with helicity 
r + is denoted as Fr_ T+ (s, cos #\y) 1001! where A (= ±i) denotes the helicity of the e~, s is the square 
of the centre-of-mass energy and cos #\y is the W production angle in the centre-of-mass frame. The 
SDM (p) elements are normalised products of these amplitudes given by Q : 

(s, cos W ) = ^ (2) 



J T—T—TJ r TJ r 



Z^At + T_ \ r T-T. 



The normalisation ensures that ^ T T p T _ T _ T+T+ = 1. The p matrix is Hermitian, giving 80 indepen- 
dent real coefficients which may be experimentally measured. The diagonal elements are defined as 
the subset of elements where r_ = tL and r + = t+. These diagonal elements are purely real and are 
equivalent to the probability of producing a final W + W~ state with helicities r + and r_ respectively. 
The off-diagonal elements represent interference terms between the different helicity states of each W. 

In the context of a limited sample of events it may not be possible to measure all of the components 
independently. It is then useful to consider the single particle density matrix for either the W~ or W + , 
formed by summation over the helicity states of the other W. For example the W~ matrix elements 
p^_ T i_ are obtained by summation over r + = r+: 

PV-TL ( S i COS #W) = ^2 Pt-tLt + t+ (s, cos 0w) (3) 

The Hermitian matrix p has eight independent real coefficients, and satisfies the normalisation con- 
straint ^2 T _ p^_ T _ = 1- The three real diagonal elements represent the relative production probabilities 
for a final state W~ with a particular helicity. 

The constraints of CPT and CP-invariance impose additional symmetries on the density matrix 
at tree level [g]. CPT-invariance imposes the conditions: 

Re«; 2 ) - He(p^_J=0 (4) 

Im«; 2 ) + Im(p H ;;_ T2 )=0 (5) 

CP invariance imposes the condition: 

Im(p^) - lm(p^_ r2 )=0 (6) 

Thus CPT and CP-conservation together dictate that all coefficients from p^ T2 are real. Deviations 
from the validity of equation (y) would represent an unambiguous signal for CP-violation at tree level. 



Deviations from the equality of equation (||) represent a signal of higher-order (loop) effects beyond 
tree level §. 

The exclusive differential cross-sections for the production of a W with transverse (T) or longitu- 
dinal (L) helicity are obtained by weighting the total differential cross-section by the relevant elements 
of the single particle p matrix: 

dux dcr 

(/>++ + P- 



dcos #w dcos #w 

and the corresponding total cross-sections are given by: 

" +1 da 



dcos #w dcos #\y 

da L da 

(Poo) -r— ~ a (7) 



O-T = / (p+++p— ] 



dcos #w 



dcos #w 



.! (poo) ta% dcos% - (8) 

The single W SDM elements can be obtained from measurements of the properties of the W decay 
products by the application of suitable projection operators K TT > |l],|2| which assume the V— A coupling 
of the W to fermions. The W decays are characterised by the polar and azimuthal angles of the decay 
fermion in the W rest frame, 9* and <f>* respectively. In this analysis we consider only the semi-leptonic 
event class, WW — ► qqiv, where one W decays to a lepton (e, \x or r) and a neutrino and the other to 
two hadronic jets. In this case the values of 9* and <f>* of the lepton may be determined unambiguously 
and the corresponding single W SDM elements are given by: 



da(e + e- -» W+W" 



' w- 

Ptt' 



dcos #w 

1 f dcr(e+e- -»■ W + £~P) 



Br(W- ->£-v) J dcos (9 w dcos 9*dcf>* 



A TT ,(9*,(P*)dcos9*d(f)* (9) 



In section || we describe how this expression is realised as a sum over the events observed in the data 
sample. 

In the case of the hadronically decaying W differentiation between the particle and anti-particle 
decay products is very difficult. However, certain projection operators are symmetric under the trans- 
formation cos 9* — ► — cos 9*, (jf — > (ft* + 7T, i.e. the interchange of one of the decay particles from the W 
boson with the other decay particle, for example the u-type quark with the d-type in a hadronically 
decaying W boson. This means that a number of combinations of the SDM elements may still be 

extracted from W — ► qq decays: these include p ++ + p and poo- Both W bosons in the event can 

therefore be used to measure the polarised cross-section, since only these terms appear in equations 
©■ 

Returning to the full two particle SDM, the analogous expressions to (0) but describing the dif- 
ferential cross-section for both W bosons in the pair being transversely polarised (TT), both being 
longitudinally polarised (LL), or one of each polarisation (TL) are: 

do"TT / x dcr 

(P++++ + P++— + P—++ + P- 



dcos#\y '' dcos6*w 

d<7 LL da 

(Poooo)^ -r- (10) 



dcos 0\y dcos #w 

dcr T L , . dcr 

IP++00 + P 00 + Poo++ + Poo — )- 



dcos #\y dcos ^w 

6 



and the analogous expression to @ is 

ckr(e+e- -► W+W~ 



dcos^w T T - T+T + 

1 f d<r(e + e~ — ► qqlv) 



Br(W -> qq)Br(W -> Iv) J dcos 6» w dcos 0* d^ldcos 6* + d<f)* + 
x A t _ tL (9*_ ,cj)^)A T+T , + {9* + ,(j)* + )dcose*_d<p*_dcos6* + d(l)* + 
where the integral is now over the decay angles of both W bosons. 

3 The Data Sample and Monte Carlo Simulated Events 
3.1 The Data Sample 



(11) 



The W-pair data used in this analysis were collected by the OPAL [18[| detector at LEP. The accepted 
integrated luminosity in 1998, evaluated using small angle Bhabha scattering events observed in the 
silicon tungsten forward calorimeter fl9|| , is 183.14 ± 0.55 pb _1 pQ| . The luminosity-weighted mean 
centre-of-mass energy for the data sample is yfs = 188.64 ± 0.04 GeV. 

In this analysis we use only the W-pair events decaying to the qqeu, qqjiv and qqrv channels. These 
qqlv events were first selected using the likelihood selection described in |||, E(J,El|]. The selection is 
designed to optimise the rejection of Z°/7 — > qq and four-fermion backgrounds. A total of 1252 qqlv 
candidates was selected at this stage. Monte Carlo studies show this selection is about 88% efficient. 

Further cuts are applied in order to obtain a sample of well reconstructed events. A full description 
can be found in [|10[ . A brief overview of the procedure is given here. For the qqef and qq^iu events 
a well reconstructed lepton track is required. For the qqrf events, either one track or a narrow jet 
consisting of three tracks is assigned as the r decay product. 

Kinematic fits are now applied to the events. For the qqef and qq^iv events a one-constraint fit is 
applied that requires energy-momentum conservation and neglects any initial state radiation. Events 
are accepted if the fit converges with a probability larger than 0.001. An improvement is made to 
the resolution of the angular observables by performing a second kinematic fit which constrains each 



reconstructed W mass to the average value measured at the Tevatron [22|. Events for which the fit 
converges with a probability of at least 0.001 are accepted; for other events we revert to the previous 
fit. 

For the qqri^ events a different kinematic fit is applied f9|; here the reconstructed mass of the 
two W bosons is required to be equal. In order to obtain a W mass from the leptonic part of these 
events, it is necessary to assume that the direction of the visible part of the r decay approximates the 
direction of the r lepton. The fit is required to converge with a probability of at least 0.025. This 
cut rejects 41% of the background. It also rejects about 14% of the signal events, but preferentially 
rejects events where the r decay products are not correctly identified. 

After application of all selection cuts 1065 qqiu candidates remain, with 359 in the qqef channel, 
386 in the qq^iv channel and 320 in the qqrf channel. The W production and decay angles for these 
events are shown in figure |l[ The background contributions to the data sample are estimated using 
Monte Carlo samples. The resulting contaminations are expected to be four-fermion (non-CC03) 3.1%, 



Z°/7 — ► qq 2.1%, two-photon (multiperipheral) 0.2% and CC03 four-fermion (WW — > qqqq and tvlv) 
0.5%. 



3.2 Monte Carlo 

A number of Monte Carlo models are used to provide estimates of efficiencies and backgrounds as 
well as the expected W-pair production and decay angular distributions for different TGC values. 
The Monte Carlo samples are generated at y/s = 189 GeV. All Monte Carlo samples mentioned 
below are passed through the full OPAL simulation program |23] and then subjected to the same 
selection and reconstruction procedure as the data. The ERATO [|13J| Monte Carlo program is used in 
this analysis because it can generate samples of four-fermion Monte Carlo events with non-Standard 



Model values of CP-violating anomalous couplings. The EXCALIBUR [24] Monte Carlo is also used 
extensively in this analysis for correction of detector effects and as the main tool to compare to the 
189 GeV data. A comparison between ERATO and EXCALIBUR angular distributions and spin 
density matrix elements was undertaken and no statistically significant differences were seen between 
them. The other four-fermion generator used to calculate background contributions and systematic 
uncertainties is grc4f |2J| . The WW generators used to calculate theoretical predictions for the W-pair 
polarisations and to calculate systematic uncertainties are EXCALIBUR, HERWIG |26 1 , PYTHIA p?]] 



and KORALW |2S(| . The background Z°/7 — > qq samples are generated using PYTHIA. For samples 
of two-photon (multiperipheral) background PYTHIA, PHOJET pi and HERWIG were used. 



4 Measurement of W Boson Polarisation 



4.1 Experimental Method 

Measurements of the production cross-sections for each of the transverse and longitudinal states of the 
W bosons are obtained from summations of the SDM elements obtained from the qq£f data sample, 
after correcting for detector acceptance, resolution effects and background contamination. 

To calculate the SDM elements the events are divided into eight equal bins, k, of cos#w I n each 
bin the SDM elements are obtained by a summation over events, i, weighted by the relevant SDM 
projection operators as shown in equation (^). This is effectively the realisation of equation (||) as 
a summation over observed events. Nk is the number of events in bin k. In the extraction of all 
the individual W SDM elements CPT invariance is assumed, so all leptonic (hadronic) W + and W - 
decays can be combined. 

N k 

p k TT ,(co S 8 k w ) = — £A TT ,(cos0?,#) (12) 

JVfc i=i 

Estimates of the production cross-sections are derived from the diagonal SDM elements; therefore 

the projection operators only involve the polar decay angle cos 6*. It is not required that the SDM 

elements be positive definite. Both the leptonic and hadronic W decays can be used since the SDM 

combinations occurring in equation (0) are symmetric with respect to the use of the fermion or anti- 

fermion polar decay angle. Each event is also weighted by a correction factor, /, for detector acceptance 

and resolution effects: 

Nk 
Poo 



1 ' Et^T A oo(^) (13) 
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1 Nk 1 

P++ + P-- = t^f E 77^ A ++^) + A "^*)) ( 14 ) 



W 



/ 7= t ^(^ 



where A^ or is the corrected number of events in bin k: 

N k 1 

"**-£«*>■ (15) 

The factors /&(#*) are obtained from fully simulated Monte Carlo events, and are calculated as a 
function of cos#\y and cos#*. They are defined to be the ratio of the number of reconstructed events 
to the number of generated events in each bin. The correction factors have an effect of between 5% 
and 20%. The typical resolution of the measured cos#w is found to be 0.05, which is less than 20% 
of the bin width used in correcting this distribution. The typical resolution of the measured cos 0* is 
found to be 0.07, which is 70% of the bin width used in the calculation of the correction factors. 

The distribution of uncorrected events, shown in figure [IJ, also has to be corrected for detector 
acceptance and resolution effects. This is done in a similar manner to the correction of the SDM 
elements by calculating a correction factor from Monte Carlo events in a bin-wise fashion and applying 
it to the data. The estimated background contribution is subtracted from the SDM elements and the 
dcr/dcos #\y distribution. 

The polarised differential cross-sections, d<7L/dcos#\y and dax/dcos^w! are then obtained by mul- 
tiplying the binned unpolarised differential cross-section by the relevant SDM combinations, following 
equations (0). 



4.2 Experimental Results 

Figure ^ shows the individual W boson differential cross-sections obtained from the data by this 
method. The transverse and longitudinal components obtained from the leptonically decaying W and 
the hadronically decaying W are shown separately, together with the values obtained by combining 
the two. The polarisations of the W bosons in the W-pair event are correlated. The correlation is 
estimated to be 0.07, and is taken into account in the errors on the combined cross-sections. 

The fraction of each polarisation state, obtained by integrating over cos #\y and then dividing by 
the total cross-section, is given in table [l|. Once again the correlated polarisations of the two W bosons 
are taken into account when deriving the errors. The estimation of the systematic errors shown in 
tables IT] and is described in section H. 

The joint polarised cross-sections are obtained by extracting the joint SDM elements in a similar 
way to equations (|i~2|) 

1 Nk 1 

- ^E , , fl .J ^ Ar^ lept )A^* had ) (16) 



" T+T+ " N™ ^ f k (9f ept ,et had ) 

where in this case the p and A are those corresponding to the particular SDM combinations appearing 
in equation (|10|). Note that there is now a projection operator for each W, and the correction factor 
is binned in terms of both polar decay angles. Figure |3| shows the joint differential cross-sections 
obtained by this method, these cross-sections are not constrained to be positive definite. 

The fractions of each helicity state, obtained by integrating over cos #w of the W and then dividing 
by the total cross-section, are shown in table 0. These results are highly correlated. The correlations 





0"T/Ctotal 


O'L/Ctotal 


Data 

W^ qq 
All 


0.842 ± 0.048 ± 0.023 
0.738 ± 0.045 ± 0.025 
0.790 ± 0.033 ± 0.016 


0.158 ± 0.048 ± 0.023 
0.262 ± 0.045 ± 0.025 
0.210 ± 0.033 ± 0.016 


Standard Model Expectation 

W^ qq 
All 


0.746 ± 0.006 
0.741 ± 0.006 
0.743 ± 0.004 


0.254 ± 0.006 
0.259 ± 0.006 
0.257 ± 0.004 



Table 1: The fractions of transversely and longitudinally polarised W bosons. The expected values 
are from generator level EXCALIBUR Monte Carlo, the errors being statistical only. The first error 
on the measured values is statistical and the second is the systematic uncertainty. 

between each helicity fraction, calculated from the data, are found to be: (Jtt'-&ll = 0.67±0.02, 
a TT :a TL = -0.93±0.01 and a LL :a TL = -0.89±0.01. 

The individual W polarised differential cross-sections show good agreement with the Standard 
Model predictions, as do the overall fraction of each individual W polarisation, as seen in table |j. 
The W pair polarised cross-sections show less good agreement, but OLL/ctotai is within two standard 
deviations of the Standard Model prediction and the other two are just over two standard deviations 
away. The x 2 value for the three measurements compared to the Standard Model expectations, 
including systematic uncertainties, is 4.7. This corresponds to a x 2 probability of 10%. 





Measured 


Expected 


Cr TT /<7 to tal 


0.781 ± 0.090 ± 0.033 


0.572 ± 0.010 


CLL/Ototal 


0.201 ± 0.072 ± 0.018 


0.086 ± 0.008 


OTL/Vtotal 


0.018 ± 0.147 ± 0.038 


0.342 ± 0.016 



Table 2: The fraction of each helicity combination of WW pairs. The expected values are calculated 
from Monte Carlo studies. The first errors on the measured fractions are statistical and the second 
are systematic. The errors on the expected fractions are statistical only. 



5 Measurement of Anomalous Couplings and Test of CP-Invariance 

Measurements of anomalous coupling parameters are obtained through a comparison of distributions 
of SDM elements obtained from the leptonically decaying W bosons in the qq^f data sample with 
distributions obtained from fully simulated Monte Carlo samples. In contrast to the polarised cross- 
section measurements, neither data nor Monte Carlo events are corrected for detector or acceptance 
effects. Instead the experimentally observed SDM distributions are compared directly with those for 
fully detector simulated Monte Carlo events. Monte Carlo samples for arbitrary TGC values are 
obtained by a re-weighting technique applied to a large Standard Model sample. A x 2 minimisation 
procedure is then used to find the simulated data set which best fits the data and hence obtain the 
best fit TGC parameter value. 



The events are again divided into eight equal bins of cos#w Equation (12) is used to calculate all 
six real and three imaginary SDM elements directly. The SDM distributions are shown in figure Q. 

This method is used to calculate both the CP-conserving and the CP-violating couplings. However, 
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when performing the fit for CP-conserving couplings only the six real SDM coefficients are used, and 
in the case of CP violating couplings all nine SDM coefficients are included. The effect of correlations 
between SDM elements is included. Although there is no correlation between cos#w bins, within 
each bin the measurements of the SDM elements are highly correlated. The correlations are obtained 
directly from the data. 

The Monte Carlo sample used in the fitting procedure is generated for Standard Model couplings 
using EXCALIBUR. All TGC-dependent four-fermion diagrams are included in this sample. SDM 
distributions for arbitrary couplings are obtained by re-weighting the fully simulated EXCALIBUR 
Monte Carlo events. Technically this was achieved by re-weighting events with the matrix elements 
from [l|. The five characteristic W production and decay angles are constructed from the original 
four-momenta of the primary fermions in the simulated EXCALIBUR events. This treatment neglects 
the effects of four-fermion background, but these are checked explicitly for CP-conserving couplings, 
and found to be negligible. The expected contribution from four-fermion background events produced 
in channels unaffected by triple gauge couplings, and those from other backgrounds, are added to the 
EXCALIBUR sample. 

As the SDM elements are normalised to the number of events in each cos #\y bin, information from 
the production angle of the W is effectively removed from the SDM fit. In order to include this an 
independent x 2 is derived from the comparison of the shape of the cos #\y distributions and this is 
incorporated into the fit. No information from the overall cross-section is included. 

The method of re-weighting and fitting is tested by performing fits to large samples of fully simu- 
lated four-fermion Monte Carlo samples generated with various anomalous couplings as well as others 
generated with Standard Model couplings. These bias tests show that the extracted coupling values 
are consistent in all cases with the generated values. The reliability of the statistical error is also tested 
by performing a fit to a large number of subsamples of the Monte Carlo, each with the same statistics 
as the data. It is found that in the case of the Standard Model Monte Carlo at least 68% of the fitted 
values of the couplings fall within A^ 2 =l of the Standard Model values of the couplings. Similar tests 
are performed on samples of Monte Carlo with anomalous couplings. Statistics are lower, but in all 
cases the results are consistent with the expectation, as seen for the Standard Model samples. 

For the CP-violating couplings, including statistical uncertainties and the systematic uncertainties 
described in section 0, we find: 



and for the CP-conserving couplings we find: 



K z — 


-o 2n +0 - 10 


94 = 


— n D9+ 0,32 


K = 


-0.181°;^ 


find: 




Ak 7 = 


n 99 +0.29 
u - zz -0.24 


A 5l = 


-o.o3i8:» 


A = 


-0.08+0.^ 



These latter results are consistent with those measured in [1C]. All couplings are set to their Standard 
Model values except the coupling being measured and those related to it via the SU{2)\ J x U{l)y 
symmetry. A full breakdown of the results for the CP-violating couplings is shown in table ||. 

The x 2 plots for all the fits including systematics can be seen in figure ||. The double minima in 
the x 2 curves for the CP-violating couplings derived from the cos#w distribution alone reflect that 
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the cos 0\y distribution is sensitive only to the absolute magnitude of the coupling. The same is true 
for the real SDM elements resulting in the double minima in the SDM x 2 curve for hz- It is only 
the imaginary parts that lift the degeneracy. It is evident that the CP-violating couplings are better 
constrained by the SDM elements than by the W boson production angular distribution, whereas for 
the CP-conserving couplings the converse is true. 

A further test of tree level CP violation is given by comparing the imaginary SDM elements from 
the W + to those from the W~, as in equation (y). These comparisons, which give a model independent 
test of CP-violation in the triple gauge coupling, are shown in figure ^. Deviations from zero would 
only be due to CP-violation at tree level. No deviations are seen, complementing the results of the 
measurements of the CP-violating TGCs. 

Also shown in figure are the combinations of imaginary SDM elements that test for effects beyond 
tree level, as described by equation (||). Any deviation from zero in these plots could only be caused 
by either loop effects or CPT-violation; all CP-violating tree level effects cancel out. The data are 
consistent with no effect. 



6 Systematic Uncertainties 

Systematic uncertainties in the measurements are calculated as described below. The individual 
contributions to the errors on the polarised cross-section fractions are given in table |j. For the 
measurements of the TGCs, all error sources listed below are included except that from the detector 
correction. For all TGC measurements the systematic errors are included as extra uncertainties on 
the contents of each bin of the SDM element and cos #\y distributions. Each systematic uncertainty 
is taken to be uncorrelated between bins in both cos#\y and the different SDM element distributions. 

• Jet reconstruction: Uncertainties in the modelling of jet reconstruction are estimated by varying 
the reconstruction in a Monte Carlo sample. The resolutions of the three jet parameters (energy, 
cos 8j, 4>j) are varied by 10%, and the energy is shifted by 5% to account for systematic uncertainties. 
This is done for both the quark jets and the r jets. The size of the variations for the quark jets is 
determined from extensive studies of back-to-back jets at Z° energies. A possible systematic shift in 
the reconstructed direction of the boson has been estimated using radiative Z°/7 — ► qq events. The 



possible shift in | cos #w| was found to be less than 0.01 [10|. The same uncertainties are taken on the 
r jets as for the quark jets. 

• Hadronisation: Uncertainties due to the hadronisation model are estimated by comparing Monte 



Carlo fragmented with HERWIG5.9 to Monte Carlo using the JETSET7.4 |27j hadronisation scheme. 
Both Monte Carlos have been tuned to OPAL data. Variations in the calculated polarised cross- 
sections, SDM elements and cos#\y distribution between the different samples are taken as the sys- 
tematic uncertainty. 

• Monte Carlo generators: The modelling of the four-fermion production processes is assessed by 
comparing SDM elements and cos #\y distributions from the ERATO and grc4f Monte Carlo programs 
to those from EXCALIBUR. All these generators have a different calculation of the matrix elements 
and a different treatment of ISR. The calculated values of the individual W and W-pair polarised 
cross-section fractions are also compared. Variations in the calculated polarised cross-sections, SDM 
elements and cos #\y distribution between the different samples are taken as the systematic uncertainty. 

• Background simulation: Possible systematic effects due to the simulation of the dominant Z°/7 
background are accounted for by replacing the PYTHIA Monte Carlo with a HERWIG sample. The 
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two-photon background is removed and doubled as well. Variations in the extracted polarised cross- 
sections, SDM elements and cos #w distribution are taken as the systematic uncertainty. 

• Lepton ID: The efficiency of the lepton identification in the data compared to that in the Monte 
Carlo is investigated as a function of the polar angle and energy of the lepton. The Monte Carlo is 
weighted to adjust these distributions to the data. Variations in the calculated polarised cross-sections, 
SDM elements and cos #w distribution between the samples before and after weighting are taken as 
the systematic uncertainty. 

• Detector effect correction: The use of Standard Model Monte Carlo to correct for detector 
effects could introduce some bias. This is tested for by comparing the helicity fractions calculated 
from generator level non-Standard Model Monte Carlo with those calculated from the same sample 
after full detector simulation and correction using correction factors determined from Standard Model 
Monte Carlo. This test is done with six samples of Monte Carlo each with one of the couplings Ak 7 , 
Agf, A, set to ±1. 

• Lepton energy scale and charge misassignment: Uncertainties on the lepton energy scale are 
tested by shifting the lepton energy by 0.3%. These are found to be negligible. The momentum 
resolution of electrons and muons and their charge misassignment are investigated by varying the 
resolution in Q/pt by 10%. Here Q is the lepton charge and pt is the transverse momentum with 
respect to the beam direction. 



7 Conclusion 



W-pair events with one leptonically and one hadronically decaying W are analysed to extract the 
polarisation properties of the W boson and test CP invariance. Both the individual W, and W-pair, 
polarised cross-sections have been measured. The individual W polarised cross-sections are well- 
described by the Standard Model predictions, both inclusively and separately for the leptonically and 
hadronically decaying W bosons (table ||). The results are consistent with those measured by L3 at 
the same centre-of-mass energy |15[], with the results presented here being more precise. 



The W-pair polarised cross-sections are measured for the first time at LEP. All the results are 
found to be consistent with the Standard Model expectations (table ||). 

Triple gauge couplings are extracted using information from the leptonically decaying W, together 
with the W production polar angle. They are found to be consistent with the Standard Model predic- 
tions. The CP conserving couplings have been measured elsewhere with the same data sample [|1Q|| , 
and serve as a consistency check for this analysis. The CP violating couplings k z , A z and g\ are 
measured for the first time with OPAL data, table y, and constraints are thus placed on possible new 
CP-violating contributions to the W production and decay processes. 
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Fit 


K z 


at 


A z 


SDM Elements 
cos6*w 


n i q+0.08 
-u.iy_ 007 

n i Q+o- 46 


o oo+ - 21 

u.uu_ 020 

7 +0 - 4 
u -'-1.8 


-0 12+ - 17 
-0.29t° il 


Combined 

Expected Stat. Error 


_n iq+0-06 
±0.11 


ni+ - 22 

±0.19 


-0 19+ - 18 
±0.12 


Final Fit 

Including Systematics 


-0 20+ - 10 


— D9 +0 - 32 


-0.18t° H 



Table 3: Measured values of the CP-violating TGC parameters. Both the SDM elements for the 
leptonically decaying W and the cos #\y production distribution in qqiu events from the 189 GeV data 
are used in the calculation. Errors are statistical only except in the case of the final combined fit. 
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Systematic 


TT 


LL 


TL 


T,L 


Jet Reconstruction 
WW^ qqiu 

W^ qq 


0.013 


0.004 


0.008 


0.007 
0.005 
0.011 


Hadronisatation 
WW^ qqiu 
W-»&/ 

W^ qq 


0.016 


0.002 


0.014 


0.009 
0.003 
0.021 


MC Generator 
WW^ qqiu 
W^iu 
W^ qq 


0.015 


0.013 


0.028 


0.001 
0.004 
0.006 


Background 
WW^ qqiu 
W^iu 
W^ qq 


0.004 


0.001 


0.004 


0.002 
0.003 
0.004 


Lepton id 
WW^ qqtu 
W^iu 
W^ qq 


0.017 


0.005 


0.014 


0.010 
0.017 
0.003 


Detector Effect Correction 
WW^ qqtu 

W^iu 
W^ qq 


0.004 


0.008 


0.012 


0.002 
0.009 
0.005 


Charge Misassignment 
WW^ qqtu 

W-»&/ 

W^ qq 


0.007 


0.001 


0.008 


0.005 
0.011 
0.001 


Total 
WW^ qqiu 

W^iu 
W^ qq 


0.032 


0.017 


0.038 


0.016 
0.023 
0.025 



Table 4: The contribution to the systematic uncertainty on the polarised cross-section fractions from 
different sources. 
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Figure 1: Distributions of the kinematic variables cos#w> cos 9f, cos#* 



and 



as obtained from 



the qq^ events. The points represent the data. The histograms show the expectation of the Standard 
Model and the cases Agf = ±0.5. The shaded histogram shows the non-qqii/ background. In the case 
of the W + decaying to the lepton, the value of <f>* is shifted by -k in order to overlay the W + and W - 
distributions. 
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Figure 2: The polarised differential W production cross-sections, a) is the differential cross-section for 
transversely polarised leptonically decaying W bosons in the W boson pair and b) is for longitudinally 
polarised, c) and d) are as a) and b) but are shown for the hadronically decaying W boson in the 
pair, e) and f) are the combinations of the leptonically and hadronically decaying W. Overlaid are 
the predictions for the Standard Model (solid line) and CP-violating anomalous couplings A z = —0.5 
(dotted line), k z = +0.5 (dashed line). The dotted-dashed line in b) and f) shows the zero. The errors 
include both statistical and systematic uncertainties. 
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Figure 3: The differential cross-sections for the W bosons of different helicity states, a) is the corrected 
total W production differential cross-section, b) is for production of pairs of transversely polarised W 
bosons, c) is pairs of longitudinally polarised W bosons and d) is for one of each. Overlaid are 
the predictions for the Standard Model (solid line) and CP-violating anomalous couplings A z = —0.5 
(dotted line), k z = +0.5 (dashed line). Only the solid line is visible on plot c) because this distribution 
is insensitive to the CP-violating couplings. The dotted-dashed line on plots c) and d) shows the zero. 
The errors include both statistical and systematic uncertainties. 
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Figure 4: The individual W SDM elements extracted from the leptonically decaying W boson in the 
qqiis data events at 189 GeV. The points are the OPAL data . The histograms represent the Monte 
Carlo predictions calculated from fully detector simulated Monte Carlo. The solid line shows the 
Standard Model expectation and the dotted (dashed) line that for Agf = +0.5 (—0.5). CPT invariance 
has been assumed in calculating these elements. The errors are both statistical and systematic. 
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Figure 5: The x 2 plots for the fits to the CP-conserving and CP-violating anomalous couplings. For 
the CP-conserving couplings the dashed line is the fit to just the 6 real SDM elements, for the CP- 
violating couplings it is the fit to all 9 SDM elements. The dotted line is the fit to just the cos #w 
distribution. The solid line is the combined fit. All fits include systematic uncertainties. 
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Figure 6: The combination of imaginary SDM coefficients sensitive to CP-violation at the triple gauge 
boson vertex for the 189 GeV data (three left side plots). Overlaid are the Born level predictions for 
the Standard Model (solid line) and also anomalous CP-violating couplings A z = —0.5 (dotted line) 
and k z = +0.5 (dashed line). Also shown are the combination of imaginary SDM coefficients sensitive 
only to loop effects (any deviation from zero) from the 189 GeV data (three right side plots). In all 
cases the errors include both statistical and systematic uncertainties. 
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